Assessing energy functions for flexible docking

A good docking algorithm requires an energy function that is selective, in that it clearly differentiates correctly docked structures from misdocked ones, and that is efficient, meaning that a correctly docked structure can be identified quickly. We assess the selectivity and efficiency of a broad spectrum of energy functions, derived from systematic modifications of the CHARMM param19rtoph19 energy function. In particular, we examine the effects of the dielectric constant, the solvation model, the scaling of surface charges, reduction of van der Waals repulsion, and nonbonded cutoffs. Based on an assessment of the energy functions for the docking of five different ligand]receptor complexes, we find that selective energy functions include a variety of distance-dependent dielectric models together with truncation of the ̊ nonbonded interactions at 8 A. We evaluate the docking efficiency, the mean number of docked structures per unit of time, of the more selective energy functions, using a simulated annealing molecular dynamics protocol. The largest improvements in efficiency come from a reduction of van der Waals repulsion and a reduction of surface charges. We note that the most selective potential is quite inefficient, although a hierarchical approach can be employed to take *Present address: Computational Chemistry and Molecular Structure Research, Lilly Research Laboratories, Lilly Corporate Center, Indianapolis, IN, USA Also: Department of Chemistry, University of Warsaw,